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Pulse  Detonation  Rocket-Induced 

MHD  Ejector  (PDRIME) 

-  Energy  extract  from  exhaust  flow  by 

MHD  generator 

-  Seeded  air  stream  acceleration  by 

MHD  accelerator  for  thrust 
enhancement  and  control 

Alternative  concept:  Magnetic  piston 

-  During  PDE  blowdown  process,  MHD 

extracts  energy  and  applies  to 
chamber  gas  via  a  magnetic  piston  or 
MHD  generator  to  maintain  constant 
chamber  pressure  and  temperature 


Cs  Inlet  Seed  Shock 


b) 


Magnetic  Piston  MHD  Generator 
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Simplified  Physical  Models 


Euler  fluid  transport  with  WENO 

Ut  +  F(U)X  +  G(U)y  =  S(U)  +  M(U) 
One-step  kinetics  (ionization  of  cesium) 

Cs  +  e  +  M  =  Cs+  +  M 


Fixed  electromagnetic  field  line 

-  Lorentz  force 


Fl  =  J  x  B 

J  =  o(E  +  u  x  B) 

MHD  energy  source  term 

1 2 

J  *  E  =  —  +  u  *  (J  x  B) 

G 


MHD  Generator  MHD  Accelerator 


(Nozzle)  (Bypass  tube) 
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Numerical  Simulations  of  PDRIME 


Temperature  contour  of  a  PDRIME  operating  at  25  km, 
chamber  temperature  of  3000K  and  flight  Mach  number 
of  9 

Total  impulse  per  cycle  at  30  km  and  bypass  area  per 
unit  length.  Initial  chamber  of  3000  K  and  bypass  length 
is  6  m.  The  chamber  is  seeded  with  0.5%  Cs  and  the 

bypass  tube  is  seeded  with  0.1%  Cs  (by  moles) 
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Detailed  Physical  Models 


•  Euler  fluid  transport  with  high-order  shock  capturing  schemes:  MP5, 

WENO,  etc.  in  thermochemical  nonequilibrium  (2T  model). 

•  Complex  chemical  kinetics  for  combustion  processes 

•  Collisional-Radiative  kinetics  for  electronic  excitation  and  ionization. 

•  Electromagnetic  field  coupling  via  ideal  and  resistive  MHD 


\ 

r  Pi 

■ 

■ 

\ 

'  pivn  j 

■ 

■ 

■ 

Pn 

■ 

Pn^n 

pu 

puvn  +  nxp0  -  BxBn/p0 

pv 

pvvn  4-  riyPo  ByBn/(.i0  : 

pUf 

.  Fn 

pwi'n  +  n , p0  -  BzBn//.l0 

^n^x  ^Bn 

By 

VfiSy  VBji 

Bz 

v nBz  wBn 

E 

vn(E  +  p0)  -  Bn(v  -  B)/jj,0 

Ph  j 

L  V<nSe  J 

A 


Distribution  A:  Approved  for  public  release;  Distribution  unlimited 


HPC  Capability:  Graphic  Processing  Units 


Graphic  Processing  Units 

-  Graphic  processing  units  containing  a  massive 

amount  of  processing  cores  designed  for  graphic 
rendering 

-  GPU  is  faster  than  CPU  on  SIMD  execution  model. 
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Shared  Memory 

-  GPGPU  programing  is  very  identical  to  traditional 

C/C++  (or  Fortran).  Object-Oriented  features  are 
available  in  CUDA. 

-  GPU  is  much  cheaper  than  CPU.  The 

performance/dollar  ratio  for  GPU  is  at  least  16 
times  higher  than  CPU. 
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Hoffman  2  ARK  gpu  queue  configuration 

•  6  computing  nodes  (3  Tesla  Fermi  M2070  GPUs/node) 

•  Supports  MPI,  CUDA,  OpenCL,  etc. 
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High-order  Simulations  of  Fluid  Flows 


•  Standard  test  cases  for  ideal  gas  flow: 


Forward  step  problem  (800,000  cells) 


Rayleigh-Taylor  instabilities 
problem  (1 .6  Mcells) 
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•  Chemical  kinetics 
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•  Collisional-Radiative  (CR)  Kinetics 


Induction  delay  time  for  Air  mixture  at  various 

pressure  and  temperature 
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•  Implicit  kinetics  solver 


CR  kinetics  of  electronic  states  of  atomic 
hydrogen  during  isothermal  heating  process 
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dea  MHD  test  cases 


Rotor  problem 
(16,000  cells) 


Orszag-Tang  problem 

(16,000  cells) 


Resistive  MHD 
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MHD  cloud  shock 
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Comparison  of  heating  rates  (pressure  rise)  for 
explicit  and  implicit  magnetic  diffusion  schemes 
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Flame-Shock  Coupling:  1-D  Dynamics  of 

Instabilities 


•  Oscillatory  patterns  and  distinct  instability  modes  of  a  1-D  detonation  wave  at  the  CJ 

limit  can  be  observed  via  high-order  numerical  simulation. 

-  High  frequency  mode  appears  first  and  marks  the  transition  from  a  stable  CJ  detonation 

-  High  amplitude  (low  frequency)  mode  appears  later  which  directly  couples  the  flame  speed  with  the  shock 


•  The  dynamics  of  a  detonation 

wave  can  be  explained  by  the 
strong  coupling  between  fluid 
convection  and  kinetics. 

-  Simple  model:  entropy  and  acoustic 

waves  cycles  in  flame-shock  coupling 
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•  Influence  of  traverse  waves 

-  Unstable  shock  front  due  to  interaction  with 

the  traverse  waves. 

-  The  intersection  of  the  lead  shock  and  the 

traverse  waves  produces  two  triple  points 
connected  by  a  Mach  stem  that  grows  in 
width 


Trinle— Point 


•  Cellular  pattern  for  2D  detonation 

simulation  with  complex  kinetics 

-  Characterize  a  network  of  triple  points  that 

are  connected  to  smooth  shock  fronts 
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GPU  parallelization 


Approaches 

-  Cell-based  (1  cell  per  thread)  parallelization:  EOS,  time  marching,  chemistry,  etc. 

-  Face-based  (1  face  per  thread)  parallelization:  Reconstruction,  flux,  etc. 

Optimization  strategies 

-  Coalesce  memory  access  for  optimal  DRAM  memory  bandwidth 

-  Utilize  on-chip  shared  memory  whenever  possible 

-  Reduce  memory  transfers  between  cpu  and  gpu 

-  Adapt  instruction  level  parallelism  (ILP)  with  in  GPU  kernels  by  reducing  block  occupancy  and  utilizing 

more  registers 

Message  passing  interface  (MPI)  for  multiGPU  computation 

-  Domain  decomposition 

-  Boundary  exchange 
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GPU  vs.  CPU 


Comparison  between  NVIDIA  Tesla  M2070  and  Intel  Xeon  X5650  (single  core 

comparison) 


Ideal  Gas  simulation  (no  chemical  reaction) 

-  5th  order  accurate  solutions  (MP5  and  ADERWENO) 

-  Speedup  factor  up  to  60 


Reacting  flow  simulation  (detonation) 

-  Comparison  of  two  different  mechanism:  H2-Air  and 

CH4-Air  detonation 

-  Speedup  factor  up  to  40 
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